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^ ■ Abstract 

We consider the image denoising problem using total variation regularization. This problem 
is computationally challenging to solve due to the non-diffcrcntiality and non-linearity of the 
regularization term. We propose a new alternating direction augmented Lagrangian method, 
f**- ' involving subproblems that can be solved efficiently and exactly. We compare our method with 

the split Bregman method |10j and demonstrate the superiority of our method in image restora- 
tion quality on a set of standard test images. 
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1 Introduction 

In signal processing, total variation (TV) regularization is a very popular and effective approach for 
^ noise reduction and has a wide array of applications in digital imaging. The total variation is the 

integral of the absolute gradient of the signal. Using TV regularization to remove noise from signals 
was originally proposed in [21] and is based on the observation that noisy signals have high total 
variation. By reducing the total variation of a noisy signal while keeping the resulting signal close 
to the original one removes noise while preserving important details such as sharp edges. Other 
existing denoising techniques include median filtering and Tikhonov-like regularization, H^Ht/a - := 
^j(V x u)? + (y y u)f. It is known that they tend to smooth away important texture details along 
with the noise |231 [26] . 

For a 2-D signal u £ M nxm , such as an image, the total variation [|u[|rv [21] of u can be defined 

■ anisotropically or isotropically: 

| \V x u\ + | V a u| , (Anisotropic); 
N|TV l E.m/(V^1 2 + (V^, (Isotropic). (1) 



Concisely, IHItv can be expressed as Yli=i \\Diu\\ p , where DiU 6 M 2 denotes the discrete gradient 
of u at pixel i. Hence, ||w||rv is isotropic when p = 2 and is anisotropic when p = 1. TV denoising 
(also called ROF denoising) corresponds to solving the following optimization problem, 

min AIIuIItv + r||« - b\\ 2 , (2) 
u 2 
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where b E ^™ xm j s th e noisy image, and the solution u is the denoised image. || • || without a 
subscript denotes the ^-norm. Let us first focus on the case of anisotropic TV regularization. 
The isotropic TV model will be considered in Section 12.31 We assume that all 2-D images are in 
column-major vectorized form. Assuming that the equivalent one-dimensional index of is k 
and 1 < i < n, 1 < j < m, the elements of Vii are given by 

[V«]v = ( Uk+1 ~ Uk ) = ( Z XU V (3) 

1.1 Related Work 

Due to the non-differentiability and non-linearity of the TV term in problem ([2]), this problem is 
computationally challenging to solve despite its simple form. Hence, much effort has been devoted 
to devise an efficient algorithm for solving it. A number of references are provided in Section 1 of 
|10j . In addition, Chambolle's Algorithm [14] solves problem ([2]) with the isotropic TV-norm. In |26|, 
variable-splitting and penalty concepts are applied to solving the TV-based deblurring problem 

mm\y^\\Diuh+ -\\Ku-bf. (4) 
u £ — ' 2 

i 

Specifically, an auxiliary variable w% 6 M 2 is introduced for each pixel to decouple Diu from the 
non-differentiable term Y2i II " 111; yielding an approximation to problem Q 

minA\^ \\wi\\i + t; \\ w i ~ D M\ 2 + \\\ Ku ~ b\? '■ ( 5 ) 

w,u * — • Z * — * A 

i i 

Problem ([5]) is then minimized alternatingly with respect to w and u, with a continuation scheme 
driving the penalty parameter f3 gradually to a sufficiently large number. This method is extended in 
|281 130] to solve the multi-channel (color) image deblurring problem. In [30], the TV regularization 
with 1-norm fidelity (TVL1) model 

minA> II-Dj^II + \\Ku — b\U 

i 

is considered. The same approach has also been applied to reconstruct signals from partial Fourier 
data in the compressed sensing context [31] . A downside to this quadratic penalty approach is that 
when P is very large, problem ([5]) becomes ill-conditioned and numerical stability becomes an issue. 

Our algorithm is most closely related to the split Bregman method [10] , which is an application 
of variable splitting to the Bregman iterative regularization method [16J. The Bregman iterative 
regularization method was first introduced in [16] as a better alternative (iterative) approach to the 
TV denoising problem than directly solving problem ([2]). The Bregman distance associated with a 
convex function E(-) between u and v is defined as 

D p E (u, v) := E(u) - E(v) - p T (u - v), 

where p S dE{v) and dE(v) denotes the subdifferential of E(-) at the point v. The Bregman 
iteration for the unconstrained minimization problem 

minE(u) + XH(u), 
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where both functions E(-) and H(-) are convex, is 

u (fc+i) = m i n D p Ju,uW) + \H(u) 

u 

= mmE(u) - {u-u {k) ) T P [k) + XH(u), (6) 



p (k+l) = p (k) 



VH( U ( k+1) ). (7) 



If we linearize the H(u) term in ([6]) and add an ^-proximity term ^\\u — u( fc )|| 2 , we obtain the 
linearized Bregman iteration |33|, [52] 

u {k+l ^ = mmD p E (u,u^) + XVH(u^fu+^-\\u-u^\\ 2 , (8) 
(fc+D = p ik) _ VH{u {k) ) _\ {u ik+i) _ u (k)y (9) 

With the introduction of an auxiliary variable d in the spirit of [26], the TV denoising problem ([2]) 
is equivalent to the constrained problem 

min \\\d\\i + R{u) (10) 

u,d 

s.t. d = $(u) 



2) 



where R(u) = ^\\u — b\\ , and <&(u) = ( ^ x ). Now, converting problem f)10|) into an uncon- 



strained problem (by penalizing \\d — ) 

min A||d||i + R(u) + -L||d - <5>(u)\\ 2 

u,d 2/1 

and applying the general Bregman iteration ©-(0) with E(u,d) = A[|d[|i + R(u) and H(u,d) = 
\\d — <!>(«) || 2 , we obtain after simplification the following specific Bregman iteration: 

( u (fc+l) )d (fc+l)) = min x\\d\\ 1+ R( u ) + —\\d-<S>(u)-rW\\ 2 , (11) 

u,d 2/i 

r (k+i) = r (fe) + ($(„(fe+l))_d( fe + 1 )), (12) 

with r'°) = 0. As is well known, the Bregman iterative algorithm (jlip -(|12 p is equivalent to applying 
the augmented Lagrangian method |13[ 117] to solve problem (|10p . In [10 1, an approximate solution 
to (jlip was proposed by alternatingly minimizing the right-hand-side of (jlip with respect to u and 
d once. This yields the following split Bregman, or equivalently, alternating direction augmented 
Lagrangian (ADAL) algorithm ( Algorithm [LI]). In this paper, we propose a different ADAL method 
by further splitting the vector d. 



1.2 Organization of The Paper 

The outline of the rest of the paper is as follows. We first describe our proposed alternating direc- 
tions augmented Lagrangian method in Section [2.21 We then discuss in Section [2.41 the difference 
between our algorithm and the split Bregman method. In Section [3j we test both algorithms on 
a set of standard test images and demonstrate the superiority of our method in image restoration 
quality over the split Bregman method. 
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Algorithm 1.1 SplitBregman 



1: Given «(°), dP\ and r®*. 

2: for fc = 0, 1, • • • ,K do 

3: «( fe+1 ) <- min UJ R(n) + ^||d (fc) - *(«) - r( fe )|| 2 
«- mm d \\\d\\i + j^\\d- <3?(u( fc+1 )) - r^> 
5- r (fe+l) ^_ r (fc) + ($( lt ( fc + 1 )) _ rfC^ 1 )) 

6: end for 

7: return vS K ^ 



2 An Alternating Directions Method 

Our strategy is to reformulate problem ((2j) into an equivalent constrained optimization problem 
and use the Alternating Directions Augmented Lagrangian (ADAL) method to solve the resulting 
problem. 

2.1 ADAL 

The ADAL method is also known as the alternating direction method of multipliers (ADMM) and 
was first proposed in the 1970s 0[9]. It belongs to the family of the classical augmented Lagrangian 
(AL) method |17 } I20 | H5]. which iteratively solves the linearly constrained problem 

min Fix) (13) 

X 

s.t. Ax = b. 

The augmented Lagrangian of problem ([13]) is C(x,j) = F(x)+j T (b— Ax) + i^\\Ax — 6|| 2 , where 7 is 
the Lagrange multiplier and \x is the penalty parameter for the quadratic infeasibility term. The AL 
method minimizes C(x, 7) followed by an update to 7 in each iteration as stated in the following 
algorithm. We denote by K the user-defined maximum number of iterations or the number of 
iterations required to satisfy the termination criteria. 

Algorithm 2.1 AL (Augmented Lagrangian method) 

1: Choose 7^. 

2: for A; = 0,1, ••■ ,K do 

3. x( fc+1 ) •(— argmin^ £(x,7^)) 

5: end for 

6: return x^ K * 



For a structured unconstrained problem 

mm F(x) = f(x) + g(Ax), (14) 

X 

where both functions /(•) and <?(•) are convex, we can decouple the two functions by introducing 
an auxiliary variable y and transform problem (|14p into an equivalent linearly constrained problem 

min f(x) + g(y) (15) 
x,y 

s.t. Ax = y. 
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Now applying Algorithm 12.11 to the above problem, the augmented Lagrangian becomes 

£>{x,V,l) = f(x) + g(y) + j T (y - Ax) + — \\Ax - y\\ 2 . 

Exact joint minimization of C{x, y, 7) with respect to both x and y is usually difficult. In practice, an 
inexact version of the AL method (IAL) is often used, where C(x,y,j) is minimized approximately 
as opposed to Line [3] in Algorithm 12.11 Convergence is still guaranteed in this case, as long as the 
subproblems are solved with an increasing accuracy |20| . 

ADAL (Algorithm 12. 2p is an extreme case of IAL in that it finds the approximate minimizer of 
C(x, y, 7) by alternatingly optimizing with respect to x and y once. This is often desirable because 
joint minimization of C(x,y, r y) even approximately could be hard. 

Algorithm 2.2 ADAL (ADMM) 

1: Choose 7^. 

2: for k = 0,1,- •■ ,K do 

3: x( fc+1 ) <— aigmin x £(x,y^ k \^ ( - k " > ) 

4. y(k+i) ar g mm ^ C(x^ k+1 \ y, 7W) 

5. 7 (fc+i) ^_ 7 (fc) _ 1 _ y (fc+l)) 

6: end for 

7: return j/W 



Inexact versions of ADAL (IADAL), where the subproblems are solved approximately, were 
proposed in (6j [12]. The convergence of ADAL (and IADAL, with conditions on subproblem 
accuracies) has been established for the case of two-way splitting as above [6], [12] . It is also known 
that u does not have to decrease to a very small value (or can simply stay constant) in order for 
the method to converge to the optimal solution of problem (|15p [151 [3] . 

The versatility and simple form of ADAL have attracted much attention from a wide array of 
research fields. [27] applies ADAL to solve semidefinite programming problems. In signal process- 
ing/reconstruction, ADAL is applied to sparse and low-rank recovery, where nuclear norm minimiza- 
tion is involved [141 134] . and the /i-regularized problems in compressed sensing |29j . ADAL-based 
algorithms (SALSA and C-SALSA) j2[ [Q have also been proposed to solve a number of image pro- 
cessing tasks, such as image impainting and deblurring. In machine learning, ADAL and IAL-based 
methods have been successfully applied to the structured-sparsity estimation problems [18] . 

2.2 Application to TV Denoising 

We consider the anisotropic TV denoising model ([2]). As in [10J, we introduce auxiliary variables 
d x and d y for the discretized gradient components V x u and V y u respectively. Under the Neu- 
mann boundary condition, V x u = Du, where the discretization matrix D 6 I" mx " m is an upper 
bidiagonal matrix with l's on its diagonal and -l's on its super-diagonal. Similarly, \/ y u = Dv 
and v = Pu, where P is a permutation matrix so that v is the row-major vectorized form of the 
2-D image. (Recall that u is in the column-major form.) Hence, problem ([2]) is equivalent to the 
following constrained problem 

min K¥x\\l + \\dy\\x) + -\\u-bf (16) 

d,u,v Z 

s.t. d x = Du, 
d y = Dv, 
v = Pu. 
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The augmented Lagrangian of problem (I16p is 
C(d, u, v, fi) = - \\u - b\\ 2 + A(||d x .||i + \\d y \\i) + j^(Du - d x ) + j^(Dv - d y ) + ^(Pu - v) 

1 - '|2 , II n„. j \\2\ , 1 n D „. „.ll2 



+ — {\\Du - d x f + \\Dv - d y \\ z ) + —\\Pu - vf. (17) 
2^1 l\X2 



To minimize C over d, we solve 



min A(|KHi + KHi) + ll(Du - d x ) + tf(Dv - d y ) + — {\\Du - d x f + \\Dv - d y f). (18) 

d Zfl\ 

Problem (|18p is separable with respect to d x and d y , so the unique minimizer can be computed 
through two soft-thresholding operations 

d* x = T{Du + /ii7x,A/ii), 
d* y = T(Dv + /Ui7y 5 Vi), 

where the soft-thresholding operator T is defined componentwise as 

T(x, X)i := max{|xj| — A, 0}sign(xi). 

To minimize C over w, we solve 

11 1 

min-||u - b\\ 2 + 7^-Du H \\Du - d x \\ 2 + jJPu H \\Pu - v\\ 2 , (19) 

u 2 2// ; 2// L ) 

which simplifies to the linear system 

D T D + + Mi) u = fixb + D T (d x - fi llx ) + P T (^-v - /i l7 ^ (20) 

We observe that D T D is a positive semi-definite tridiagonal matrix. Since /ii and [12 are both 
positive scalars, the matrix on the left-hand-side of the above system is positive definite tridiagonal. 
Linear systems of this special structure can be solved efficiently by the Thomas algorithm in 8nm 
flops [IT] . We denote the solution to the above tridiagonal system by u(d x , v , 7^, 7 Z ). 
Similarly, the sub-problem with respect to v simplifies to the tridiagonal system 



D T D + — / ) V = D T {dy - fXJy) + ^ Z + ^~PU. 

^2 / /^2 



(21) 



Its solution is denoted by v(d y ,v,~/ y ,j z ). 

With all the ingredients of the algorithm explained, we formally state the ADAL algorithm 
below. 

2.3 The Isotropic Case 

The isotropic TV denoising model differs from the anisotropic model in the definition of the TV 
norm. In this case, we define ||u|| j^P := £Zt J (V x u) 2 + (V y u) 2 , and the optimization problem to 
solve is 

minA||u||^ + -||-u-6|| 2 . (22) 
u 2 
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Algorithm 2.3 ADAL (TV denoising) 



1: Given^), W W,d(0),A, 7 i 0) , 7 f, 7 f. 

2: for A; = 0,1, ■■■ ,K do 

5: <- v{d { y +l \u {k \-f ( y\^ ] ), the solution of fl2TJ 

6: n( fc+1 ) <- n(di fe+1) ,u( fc+1 ),7i fe) ,7l fe) ), the solution of ((20j) 
7: 7 ^ +1 V 7 « + ^(W+ 1 )-4 fe+1) ) 
8: 7 r i V 7 f + J r (^)-4 fe+1) ) 

9: 7 i fc+1) <- 7 i fc) + ^(W fc+1 > - ^ fc+1 >) 

10: end for 

11: return u^ K ' 



We observe that by the definition above, H^llyy = z~2i [V y u]i)||, which is the group lasso 

regularization on (V x u,V y u), with each group consisting of ([Vxitji, [V y it]i). We introduce the 
same auxiliary variables and linear constraints among them as in the previous section. As a result, 
the two subproblems with respect to u and v are the same as (|20p and (|2ip respectively. The 
subproblem with respect to d now becomes 

njnA^ \\{[d x ] h [d y ]i)\\+^(Du - d x ) + ^{Dv - d y ) + ^-{\\Du - d x f + \\Dv - d y \\ 2 ), (23) 

i 

which is a proximal problem associated with the group ^2-norm ||d||i,2 with d x = V x u,d y = V y u, 
and the groups being defined above. The solution to this subproblem is thus given by a block 

soft-thresholding operation [25j HHJ E], <S(f ^ ^ J ( ^ J + fii f ^ , Xfi\), where the block 

soft-thresholding operator is defined blockwise as 

S(x, X)i := max{||xj|| — A, 0} 1 , 



and Xi is the i-th block of x, i.e. ([Du]i, [Dv]i) in our case. 
2.4 Comparison with The Split Bregman Method 

The Bregman iterative regularization method has been shown to be equivalent to the classical aug- 
mented Lagrangian method [33]. Moreover, the alternating minimization approach for minimizing 
the subproblem (jlip makes the split Bregman method (Algorithm II. ip equivalent to Algorithm 12.21 
[24], [7j [22] applied to the constrained problem 

min A(||4.||i + \\d y \\i) + -\\u - b\\ 2 

d,u 2 

s.t. d x = V x u, 

dy = X7yU. 

It is clear now that the main difference between Algorithm 12.31 (ADAL) and the split Bregman 
method is the additional constraint v = Pu in problem (I16p . The split Bregman subproblem with 
respect to u (line [3J in Algorithm 11. ip can be simplified to the linear system 

(/xJ + {V T X V X + VjV„)) = lib + VM* - r x k) ) + V T y(d^ - r< fc >), (24) 
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Figure 1: The set of standard test images. 



whose left-hand-side matrix includes a Laplacian matrix and is strictly diagonally dominant. Solv- 
ing this linear system exactly in each iteration is expensive. Hence, one iteration of the Gauss-Seidel 
method is applied in [10] to solve f)24f) approximately. 

In contrast, the subproblems with respect to v and u ( (|2ip and (|20p ) in ADAL (Algorithm 12. 3p 
have much better structures and thus can be solved exactly in an efficient manner as we saw in 
Section 12.21 We surmised that this is a better approach for the TV denoising problem; the results 
in the next section confirmed this. 



3 Experiments 

Our ADAL algorithm (Algorithm 12. 3p was written in Matlab, whereas SplitBregman is in C with 
a Matlab interface. We ran both algorithms on a laptop with an Intel Core 2 Duo processor and 
4G memory. 

3.1 Test images 

We compared our ADAL algorithm with the split Bregman method on a set of six standard test 
images: lena, cameraman, house, blonde, mandril, and peppers. They present a range of 
challenges to image denoising algorithms, such as the reproduction of fine detail and textures, sharp 
transitions and edges, and uniform regions. Each image is 512 x 512 in grey-scale pixels and is 
denoted by uo in vectorized form. The original images are presented in Figured) 

1 Code downloaded from http://www.stanford.edu/ tagoldst/code.html. 



8 



3.2 Set-up 



We constructed the noisy image by adding Gaussian noise to the original image, i.e. b = uq + e, 
where e ~ A/"(0, a) and b is the vectorized noisy image. The quality of the denoised image in 



peak-signal-to-noise ratio (PSNR) to compare the reconstruction quality. The PSNR of an image 
u with respect to the noiseless image no, in the case where the maximum pixel magnitude is 255, 
is defined as 



Note that the PSNR is monotone decreasing with the normalized error, i.e. a higher PSNR indicates 
better reconstruction quality. We tested three cases with the anisotropic model: low noise (a = 30), 
medium noise (<r = 50), and high noise {a = 75) with a range of values for A. In order to select 
the appropriate A's for the different levels of noise, we considered both the PSNR and visual 
impression. Although the PSNR is a good indicator of reconstruction quality, it is not always 
true that a denoised image with a higher PSNR with respect to the original image is more visually 
appealing than another with a lower PSNR. Hence, for each of the three cr's, we examine in the next 
section three values of A: (1) one that gives the best visual quality, (2) one considerably smaller 
than (1), and (3) one considerably larger than (1). 

3.3 Strategy for updating /i 

We started with [i\ = 0.2 and \i2 = 1-5. We decreased [i\ in (fT7|) by a factor of (3 = 0.5 and \i2 by 
a factor of (3 = 0.1 if 77^) increased for = 3 consecutive iterations. The lower limit of ui and \ii 
was set to 10~ 4 . 

3.4 Results 

In all test cases, we ran ADAL and SplitBregman for 25 iterations, at which point both algorithms 
have converged. We report the best PSNR's obtained for the six images at each of the noise levels 
in Table [TJ The first striking observation is that ADAL outperformed SplitBregman in all the test 
cases in the table. In fact, even if we consider the other values of A, ADAL always produced a 
higher PSNR than SplitBregman. (See Figures [2j [31 and[H) We also present the convergence plots 
of the normalized errors for all test images. To avoid repetition, we show the plots only in the case 
of medium noise. Basically, the graphs show that ADAL has a very competitive convergence speed 
compared to SplitBregman. Considering that the per-iteration computational cost of ADAL is only 
slightly more expensive than the per-iteration cost of SplitBregman, the total denoising times are 
expected to be similar. 

3.4.1 A for best visual quality 

We observed in the experiments that the best visual quality for a given image and method usually 
occurred at a A slightly larger than the one giving the highest PSNR, which tended to appear 
noisier than one expects. Specifically, we found that setting A = 30 for low noise, A = 60 for 
medium noise, and A = 80 (for mandril, A = 70) for high noise gave the best visual results. The 
gaps in the normalized error and PSNR are quite small (Figure [5|), but there are subtle differences 
in the reconstruction quality as we illustrate in Figure El We have magnified the restored images 
of lena around the face. Besides a higher PSNR, the image produced by ADAL clearly appears 
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a 


Algs 


lena 


cameraman 


mandril 


blonde 


house 


peppers 


30 


ADAL 


29.299 


29.836 


25.481 


27.919 


32.094 


29.515 


SplitBregman 


29.125 


29.615 


25.115 


27.893 


31.815 


29.273 


50 


ADAL 


27.042 


27.138 


23.493 


25.735 


29.984 


27.495 


SplitBregman 


26.910 


27.062 


23.327 


25.630 


29.814 


27.342 


75 


ADAL 


25.845 


25.750 


22.157 


24.766 


28.216 


26.179 


SplitBregman 


25.731 


25.724 


22.028 


24.699 


28.215 


26.055 



Table 1: The highest PSNR's produced by ADAL and SplitBregman on each of the test images 
with different amount of noise. 



Iena.tif cameraman.tif 




Figure 2: Plots of PSNR against A for a = 30. Red: ADAL. Black: SplitBregman. 
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cameraman.tif 
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Figure 3: Plots of PSNR against A for a = 50. Red: ADAL. Black: SplitBregman. 




Figure 4: Plots of PSNR against A for a = 75. Red: ADAL. Black: SplitBregman. 
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Error plot for cameraman.tif 
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Error plot for mandril.tif 
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Error plot for blonde.tif 
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Error plot for house.tif 
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Error plot for peppers.tif 
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Figure 5: Convergence plots of normalized errors for a = 50 and A = 60. 



more appealing. In the case of mandril, there are also tangible improvements in contrast and 
texture details (e.g. the nose) in the results from ADAL. For images with more or less piecewise 
constant intensities (e.g. house and cameraman), the two algorithms rendered very similar and 
good reconstruction results. 



3.4.2 Small A 

For small A's, we set A = 10 for low/medium noise and A = 30 for high noise. Figure [7] shows that 
ADAL leads SplitBregman by a significant margin in terms of normalized error (and hence, PSNR 
as well) . We also present the final denoising results after 25 iterations on three test images in Figure 
Since small A's correspond to low regularization power, as expected the resulting images still 
appear fairly noisy. However, the images produced by ADAL clearly have better contrasts than 
those produced by SplitBregman. 



3.4.3 Large A 

In the case of large A's, we set A = 80 for low noise, A = 110 for medium noise, and A = 120 
for high noise. From the denoised images produced by both algorithms, we see that there is an 
over-smoothing effect due to the large TV regularizations. Although the normalized errors achieved 
by the two algorithms are quite close (Figure [9]) , the difference in the image restoration quality 
is apparent in Figure [10J ADAL was able to render a smooth surface in areas with a gradual 
color transition (e.g. the pepper skin and the women's faces) while maintaining sharp edges at the 
boundaries of different objects. In contrast, the images produced by SplitBregman appear patchy 
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noisy image, cj= 30 X= 30 PSNR= 18.590 ADAL, PSNR= 32.094 SplitBregman, PSNR= 31.815 




noisy image, o= 50 k= 60 PSNR= 14.165 ADAL, PSNR= 26.579 SplitBregman, PSNR= 26.435 




noisy image, a= 75 h= 70 PSNR= 10.612 ADAL, PSNR= 21.937 SplitBregman, PSNR= 21.780 




Figure 6: Comparison of reconstruction quality for best A's. The denoised images of lena have 
been zoomed in around the face. 
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Error plot for lena.lif Error plot for cameraman.tif Error plot for mandril.tif 




10 20 30 10 20 30 10 20 30 



Error plot for blonde.tif Error plot for house.tif Error plot for peppers.tif 




10 20 30 10 20 30 10 20 



Figure 7: Convergence plots of normalized errors for a = 50 and A = 10. 

and leave an impression of flatness. This is because there are unnecessary sharp edges in areas of 
color gradients, e.g. shades. 
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Figure 8: Comparison of reconstruction quality for small A's. 
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Error plot for blonde.tif 
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Figure 9: Convergence plots of normalized errors for a = 50 and A = 110. 
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noisy image, o= 30 k= 80 PSNR= 18.566 ADAL, PSNR= 24.211 SplitBregman, PSNR= 24.116 




noisy image, a= 75 X= 120 PSNR= 10.614 ADAL, PSNR= 24.679 SplitBregman, PSNR= 24.493 




Figure 10: Comparison of reconstruction quality for large A's. The denoised images of blonde 
have been zoomed in. 
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